Z 2 index theorem for Majorana zero modes in a class D topological superconductor 
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We propose a Z2 index theorem for a generic topological superconductor in class D. Introducing a 
particle-hole symmetry breaking term depending on a parameter and regarding it as a coordinate of 
an extra dimension, we define the index of the zero modes and corresponding topological invariant 
for such an extended Hamiltonian. It is shown that these are related with the number of the zero 
modes of the original Hamiltonian modulo two. 
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Topological invariants are useful tools in various fields 
in physics. In particle physics, interesting phenomena 
such as chiral and gauge anomalies, instantons, vortices, 
Skyrmions and many other nonperturbative effects are 
related to the topological invariants [HQ. In condensed 
matter physics, they also play a crucial role in the classi- 
fication of various kinds of phases of matter [§-0] ■ Some 
kinds of defects {8\ or solitons Q are also classified by 
the use of topological numbers. 

Recently, index theorems 0, [l(| for fermions coupled 
with a Higgs field in a monopole or vortex background 
fill HH have been attracting renewed interest in con- 
densed matter physics. It is due to zero-energy Majo- 
rana bound states obeying non-Abelian statistics, pio- 
neered by Read and Green [l]| in a p-wave supercon- 
ductor and by Kitaev [13, EH in quantum computations. 
Extensive studies have recently predicted them in vari- 
ous kinds of superfluids and superconductors 15-3^|. In 
a special case in which the systems have enhanced sym- 
metry (chiral symmetry), it has been shown [3fj| that the 
index theorem ensures the existence of such zero modes. 
However, for generic topological superconductors, there 
are no index theorems which relate zero modes with a 
topological invariant. 

In this paper, we propose a Z2 index theorem for Ma- 
jorana zero modes in a vortex of a generic topological 
superconductor in class D. We first introduce a particle- 
hole symmetry breaking term that depends continuously 
on a parameter. We next regard it as a coordinate of 
an extra dimension and extend the d-dimensional Hamil- 
tonian to (d + l)-dimensional one with chiral symmetry. 
This enables us to define the index of the extended Hamil- 
tonian and also the topological invariant corresponding 
to it. Since the index is equal to the number of zero 
modes of the original d-dimensional Hamiltonian mod- 
ulo two, we thus have a Z 2 index theorem for the orig- 
inal system, implying that the Majorana zero mode in 
topological superconductors in class D is also protected 
topologically. We will give concrete calculations in two 
dimensional models, but the extension to an arbitrary 
higher dimensions is straightforward. 

We investigate the Hamiltonian proposed by Fu and 
Kane [2lJ for the surface states of a topological insulator 



with a proximity effect of an s-wave superconductor, 



It = ij J dj — ji I 



(1) 



where 7 J and j a+2 (j,a =1,2) form a set of 7 matrices 
7'' in four dimensions satisfying {7^, 7"} = 2S^. We ex- 
plicitly employ 7 J = ® a 3 and 7 a+2 = 1 <g> a a , where, 
in the tensor product of the Pauli matrices, the former 
(latter) describes the spin (Nambu) space. The pairing 
potential is included as <fi = (ReA(ai), ImA(x)). Pro- 
vided that 4>i4>2li 7^ 0, the Hamiltonian ([1]) belongs to 
class D, since it has only particle- hole symmetry 

cue- 1 = -n, 

where charge conjugation operator C is defined by C = 
i"/ 2 "/ 3 K with K being complex conjugation 0, 0, l4fj| . 
This symmetry ensures that all states with nonzero ener- 
gies appear as paired states with ± energies. Let e n be an 
eigenvalue of the Hamiltonian ([I]) labeled by an integer n. 
Then, it is natural to label the paired state with the op- 
posite energy by — n, and hence, particle-hole symmetry 
can be characterized by the relation, e_„ = —e n . 
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FIG. 1: The spectrum of the Hamiltonian in the case of q — 3 
vorticity. (a) Three Majorana zero modes, marked by an oval, 
appear when fj, — 0. (b) Two of these get finite ± energies if 
a nonzero chemical potential is switched on, but an unpaired 
state is protected from it, sitting exactly at the zero energy 
due to particle-hole symmetry. 

Let us assume that there is a vortex at the origin which 
is described by <j> = A(r)(cos 0(0), sin 0(0)), where (r,8) 



2 



is the polar coordinates in two dimensions and the phase 
satisfies 6(27r) = 8(0) + 27rg with q being an integer. We 
assume that A(0~) = and A(oo) = A > 0. Then, it has 
been argued [ljj, [2(| that a Majorana zero mode appears 
for odd q, whereas for even q no zero modes are allowed. 
This can be understood from the perturbation theory 
based on a simpler model with /i = which bel ong s to 
class BDI due to additional chiral symmetry IE HQ- In 
this case of /i = 0, the exact q zero modes wave functions 
can be obtained similarly to 
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We can also apply the 
index theorem to this special model 0,0, and show 
that the index computed by the exact zero modes men- 
tioned above and the topological invariant, i.e., the wind- 
ing number of the pairing potential, indeed coincide, both 
of which are given by — q. This result matches the classi- 
fication scheme due to Teo and Kane [4l|. On the other 
hand, in the case of \x ^ 0, index theorems cannot apply 
any longer due to the absence of chiral symmetry. In- 
stead, perturbative arguments strongly suggest that even 
number of zero modes disappear in pairs with nonzero ± 
energy, whereas an unpaired state is protected to stay at 
the zero energy due to the particle-hole symmetry. These 
imply that the number of Majorana zero modes in class 
D is classified by Z2 , even or oddg 0MEII3. These 
features are illustrated in Fig. [1] 

Beyond perturbative arguments, we propose an index 
theorem valid for class D superconductors without chiral 
symmetry. To this end, we define, in the former part of 
the paper, an analytical index for a generic model ([1]) 
with a nonzero chemical potential. Firstly, we introduce 
symmetry breaking term 28, 41] to the Hamiltonian (pj 
such that 

H(t) = i^dj -fil®<7 3 + 7 a+2 a - A(r) 7b , (2) 

where 7b could be any hermitian matrix with C^C^ 1 = 
7b, which may be several products of the 7 matrices. 
One of simpler choices is 7b = 75 = (— i) 7 7 7 7^ For 
simplicity and for convenience, we will restrict our dis- 
cussions only to this case, 7b = 75. As for A(r), we 
assume that it is an odd function of r, A(— t) = — A(r), 
and hence A(0) = 0. We also assume that A(oo) = Ao is 
finite. Then, particle-hole symmetry can be generalized 
to 

ch{t)c- 1 = -n(-T). 

Provided that the spectrum of rL{r) is a smooth function 
of r, it can be labeled by the same quantum number n 
for the r = Hamiltonian such that %(r)ip n (T,x) — 
e n (T)tp n (T, x), where e n (0) = e n . From the point of view 
of the spectrum, the generalized particle-hole symmetry 
leads to the new relationship, 



-n(-r) = -e n (r). 



(3) 



In this sense, the states labeled by ±n can be still re- 
garded as the paired states. The spectral flow as a func- 
tion of t is shown in Fig. 
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FIG. 2: The spectral flow of the Hamiltonian for a generic 
case of /1 7^ 0. The n = state crosses the E — line at 
least t = 0. (a) If any other eigenvalue crosses it at a finite 
t, say, £ti(t~o) = 0, the spectral symmetry @ ensures that 
e_n(— To) — 0. (b) If the model parameters are changed, the 
above ±n modes may come not to cross the zero energy. Even 
in this case, the index changes by even integers because of the 
symmetry <(3j . 



Although the extended Hamiltonian ((5} does not have 
chiral symmetry, we can employ the spectral flow for the 
index theorem. To see this, we introduce a kinetic term 
for the parameter r and define 



n {3} = ia 



ir j dj + T a+3 ( j) a + iT 1 r 2 r 6 fi, 



(4) 



where j,a = 1,2,3. Newly defined T-matrices obey 
{r^, T^} = 2<5^. H {3) can be regarded as a Hamiltonian 
defined in three dimensions spanned by the coordinates 
xi,X2 and X3 = t. We have also introduced 03 = A 
and regarded it as a component of a generalized order 
parameter <p = (ReA,ImA,A). Note that the extended 
Hamiltonian (Q} has chiral symmetry T 7 H^ = —H^'T?, 
where T 7 = (— z) 3 r 1 • • • T 6 . Therefore, if the Hamiltonian 
(H]) has zero modes, they have definite chirality. 

The eigenvalue equation for the zero modes is rl^Q = 
0, which is given by H(t) as follows: 



■H(t)&. 



(5) 



To solve this equation, let us set $„(t, x) = 
fn(T)ip n (r,x). Then, d T $ n = (9 T /„)<£„ + f n {d T \)d\tp ni 
since the Hamiltonian H(t) depends on t only through 
A(t). Provided that d T X can be neglected in the adia- 
batic approximation and that ip n (r, x) is normalizable, it 
turns out that /„ is given by 



where / c is a constant spinor with T^f^ 1 — ±/^. If a 
given state for n ^ satisfies e„(+oo) > and e n (— 00) < 
0, which is the case of n = 1 and n = — 1 modes in Fig. [2] 
(a), the states labeled by n and —n are normalizable zero 
modes with chirality —1. Likewise, if e„(+oo) < and 
e n {— 00) > 0, ±n states are normalizable zero modes with 
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chirality +1. Now, define the index of the Hamiltonian 
H (3) by 

ind"H (3) = N + - N- (6) 

Then, these paired zero modes contribute to the index 
for Ti^ always by two. On the other hand, if e ra (±oo) 
has the same sign, such ±n states cannot be zero modes, 
since the wave functions are not normalizable. Therefore, 
these states give no contribution to the index. It thus 
turns out that n ^ modes can change the index by 
even integers. Contrary to these modes, the n = mode, 
which is always a zero mode of unless £o(oo) = 0, 
determines whether the index is even or odd. Thus, it 
turns out that the number of zero modes of H in Eq. (JXJ , 
which we denote as No(H), is given by the index of T~L^ 
in Eq. ^ as 

N (H) = mdH {3) mod 2. (7) 

So far we have discussed the index of the extended 
Hamiltonian Q and its modulo-two relationship with 
the number of zero modes of the original Hamiltonian 
([I]). The rest of the present paper is devoted to calcu- 
lations of the corresponding topological invariant. For 



the Hamiltonian with chiral symmetry such as T-l^ 3 \ it is 
possible to define the topological invariant equal to the 
index and to claim the index theorem 0, [l(| [36| such 
that 

ind U {3) = ~ J dSjJ j (x,0,o6), (8) 

where dSj is an infinitesimal surface elements at the 
boundary (x — > oo) of the Euclidean space R 3 , and J 1 is 
the axial vector current defined by 

J l (x,m,M) 

The first term in the above parentheses becomes the in- 
dex in the limit m — > when integrated over the two- 
dimensional surface at infinity. The second term has 
been introduced as a Pauli-Villars regulator in order for 
the current to be regularized [10|, |36[. It should be set 
AI — > oo after the calculations. When we calculate the 
r.h.s. of Eq. (jHJ, it may be easy to use the plane wave 
basis. Possible terms contributing to the index are 



/rl 3 k 
j-^GhiTrr (-iT j kj - iT a+3 <l> a + T^r 6 //) (K - (iT)A(K - ^T)A(K - fxT), 
(2ir) 



where 



G- l = l^Jkl + fc| + A 2 - nj +kj + A 2 

k = k) + ^ + [i 2 , r = 2i(r 2 r 6 /c 1 + r 6 r 1 /t 2 + r 1 r 2 A), a = ip'r a+3 3j0 a . 

i 



(9) 



The first step of the calculations is to take the trace 
of the r-matrices by the use of trTyr^^^ 1 ■ • • T^ 3 T U3 = 
(2i) 3 e^ 1V1 '"^ 3 ^ a . Lengthy but straightforward calcula- 
tions lead to 



J 1 = 2 3 ^ k e abc <t)adMd k( j) c 



S i3 ii 2 2 5 e ah3 \d 1 <j> a d 2 <p b 



d 3 k 
(2tt) 3 

d 3 k 
J2ny 



G K 



G 



As a next step, we carry out the integration over the mo- 
mentum kj in the above and over the space coordinates 
xi in Eq. (|8]). In particular in the latter integration, 
since the generalized order parameter (p depends cylin- 
drically on the coordinates (r, 9, r), it may be natural to 
regard the boundary of R 3 as a cylinder illustrated in 
Fig. [3] For convenience, let us divide it into two pieces 
I (r — > oo) and II (r — > ±oo). The calculations on these 



two regions are quite similar to those in [42j when \i = 0. 
However, in the case of fj, ^ 0, we have to choose ap- 
propriate Ao = A(oo) to make the calculations valid. To 
see this, let us set A — > A and A — > Ao, respectively, in 
the regions I and II. Then, we see that in the region I, 
G in Eq. ([9]) is always finite because of finite Ao- On 
the other hand, it becomes singular in the region II if 
| Ao | < |m|- Indeed, in this case, G~ x in Eq. ^ van- 
ishes at some momentum at the core of the vortex r = 0, 
where A(0) = 0. Therefore, we restrict our discussions 
to the case of |Ao| > M- After some straightforward 
calculations, we arrive at 



ind H< 3 > = ^L(|A -MMAo+Ml)^y 
= -sgn(A )g, 



d9e 



(10) 
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FIG. 3: The cylindrical surface of the integration in Eq. ((8]). 
"I" denotes the side of the cylinder r — > oo, whereas "II" 
denote two discs at r — > ±oo. 



where <j) a = </> a /Ao, the subscript a, & are restricted to 
1,2, and the ^-integration is over S 1 at r — > oo in the 
region II, which becomes q. It turns out that the index 
is basically determined only by the vorticity q, and is the 
same as the the index of ri when ji = 0, which has been 
calculated in ref. [36| as ind % = —q: The artificially in- 
troduced parameter A just change the sign of the index. 
It follows that the index is topological: It is indeed pro- 
tected from infinitesimal changes of the three parameters 
fx, q, and Ao- Eqs. (fTUj) and j7| are Z 2 index theorem for 
a generic model in class D. 

As we have mentioned, the topological invariant (|10[) is 
invalid if |Ao < \(J,\. In this case, the zero mode equation 
([5]) may not have normalizable solutions. This can be 
checked by the use of perturbations to this equation. Let 
us consider two possibilities of small parameters, [i and 
A. If one regards /i as a small parameter and applies 
the first order perturbation to /i = solutions, one has 
indeed one normalizable solution. On the other hand, 
in the case of small A, one can obtain the first order 
perturbative energy eigenvalue, but the first order wave 
function cannot be normalizable. This implies that in 
the case |A| <C the relationship between the spectral 
flow of H(t) and the zero mode of ri^ does not hold. 
Therefore, we conclude, from these observations, together 
with the calculations of the topological invariant, that 
the present formulation of the Z 2 index theorem requires 
I A 1 > 
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